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1.  Introduction 


The  presence  of  a  type  of  dynamical  behavior,  termed  chaotic,  has  been  ob¬ 
served  experimentally,  analytically,  and  numerically  in  many  diverse  non¬ 
linear  systems  [1,2].  Researchers  have  moved  beyond  the  point  of  catego¬ 
rizing  this  behavior  as  unwanted,  although  at  first  observation  it  may 
appear  to  be  noisy  and  unpredictable.  Since  the  fundamental  assertion  by 
Ott,  Grebogi,  and  Yorke  (OGY)  that  chaotic  dynamics  can  be  controlled  by 
the  use  of  small  perturbations  [3],  a  great  deal  of  research  has  been  con¬ 
ducted  that  not  only  experimentally  verifies  the  assertion  but  also  pro¬ 
poses  useful  applications  of  these  ideas.  Much  of  the  research  has  been 
centered  aroimd  controlling  and  stabilizing  the  unstable  periodic  orbits  [4] 
that  lie  within  chaotic  orbits.  Recent  theoretical  and  experimental  work 
shows  that  the  symbolic  dynamics  [5]  of  a  chaotic  system  can  be  controlled 
as  well  [6,7]. 

In  many  driven  nonlinear  systems,  changing  (e.g.,  increasing)  a  drive  pa¬ 
rameter  can  be  an  impetus  to  chaotic  behavior.  For  example,  a  modified 
Van  der  Pol  equation  used  by  Ueda  to  model  a  negative  resistor  oscillator 
[1]  has  the  form  x  +  {x^-l)x  +  x^  =  B  cos  cot.  This  oscillator  experiences 
steady-state,  periodic,  quasi-periodic,  and  chaotic  motion  depending  upon 
the  parameters  B  and  O).  Many  experimenters  have  observed  chaotic  be¬ 
havior  emerge  as  drive  parameters  change  [8-10];  those  regions  of  behav¬ 
ior  have  usually  been  considered  imdesirable.  With  control-of-chaos  tech¬ 
niques,  however,  it  is  possible  to  stabilize  any  of  the  dense  set  of  periodic 
orbits  embedded  in  regions  of  chaotic  behavior,  thus  creating  the  possibil¬ 
ity  of  more  robust  operation  in  the  chaotic  behavior  region. 

Present  control-of-chaos  techniques  based  upon  the  OGY  assertion  must 
wait  for  the  ergodic  wanderings  of  the  chaotic  trajectories  to  come  near  a 
desired  periodic  point.  This  is  necessary  so  that  (a)  there  is  a  linearizable 
region  near  the  periodic  point  for  which  to  calculate  control  perturbations 
and  (b)  the  control  perturbations  are  not  so  large  as  to  significantly  deform 
the  topology  of  the  state-space  attractor.  A  state-space  attractor  is  an  object 
in  the  state  space  to  which  the  trajectories  tend  [4]. 

For  some  systems  that  may  be  operating  near  crisis  [4,11],  it  may  not  be  de¬ 
sirable  to  wait  for  the  trajectories  to  wander  near  a  desired  periodic  point 
(control  point).  In  this  case,  it  is  necessary  to  use  the  inherent  sensitivity  to 
small  changes  to  guide  the  trajectories  near  a  control  point  as  rapidly  as 
possible.  This  is  the  notion  of  targeting  first  formally  introduced  in  a  cha¬ 
otic  system  by  Shinbrot  et  al  [12].  In  this  report,  we  discuss  a  particular 
method  of  targeting  using  small  perturbation  control  of  symbolic 
dynamics. 


5 


2.  Discrete-Time  Maps 

A  large  class  of  continuous-time  nonlinear  systems  that  behave  chaotically 
can  be  described  mathematically  by  a  three-dimensional  system  of  nonlin¬ 
ear  ordinary  differential  equations  [1].  The  solutions  to  these  systems  form 
a  continuous-time  flow  in  three-dimensional  space.  This  space  has  been 
termed  a  state  space,  simply  because  at  any  time  t  =  t^  the  solutions  x(tj, 
y(tg),  z(fg)  describe  a  point  (state  point)  in  the  space  that  also  describes  the 
state  of  the  system.  This  continuous-time  system  can  be  reduced  to  a 
discrete-time  map  by  judicious  placement  of  a  surface  that  cuts  through 
the  flow  of  state-space  trajectories.  The  points  that  are  produced  as  the  tra¬ 
jectories  pass  through  the  surface  create  a  two-dimensional  discrete-time 
mapping  of  state  points  called  the  Poincare  surface  of  section  [4].  Figure  1 
shows  a  continuous-time  flow  discretized  by  a  Poincare  surface  of  section 
placed  at  constant  z,  thus  producing  a  two-dimensional  map  described  by 
[Xn+y  y„+i]  =  y„)-  Poincare  surfaces  of  section  for  a  large  class  of 

highly  dissipative  nonlinear  systems  can  produce  approximately  one¬ 
dimensional  discrete-time  surfaces  [13,14].  If  we  focus  on  one  of  the 
discrete-time  state  points,  say  then  we  can  form  a  mapping  based  upon 
its  sequential  surface  crossings.  This  map  is  commonly  called  a  return  map, 
which  has  the  form  =  /(x„),  and  will  also  be  approximately  one¬ 

dimensional  for  a  dissipative  system.  We  focus  on  this  particular  case  and 
show  how  the  formalism  can  be  generalized  for  higher  dimensional  maps. 

A  paradigm  for  the  study  of  one-dimensional  chaotic  maps  is  the  logistic 
map.  It  is  a  single-maximum,  noninvertible  map  described  by  the  expres¬ 
sion  =  pz„(l  -  x^).  We  use  the  following  noise-perturbed  map: 

^n+\  ~  P^n 

where  p  is  a  parameter  that  governs  the  global  behavior  of  the  map,  and 
is  a  small  amount  of  uniformly  distributed,  zero-mean  noise  with  a  peak- 
to-peak  amplitude  of  10“^^  that  we  added  to  avoid  the  round-off  induced 


Figure  1.  A  three- 
dimensional  chaotic 
attractor  demon¬ 
strating  how  state- 
space  trajectories 
pierce  Poincare 
surface  of  section. 
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periodicity  inherent  in  computer  simulations  [15].  Figure  2  shows  the  map 
created  by  successive  iteration  of  equation  (1)  for  p  =  3.97.  The  points  that 
are  generated  by  successive  iterations  of  the  map  can  be  referred  to  as  the 
attracting  set. 

We  can  find  the  boimdaries  of  the  map  by  first  finding  the  state  point 
that  maps  to  the  maximum  edge  of  the  attractor  as  a  function  of  p.  At  the 
peak,  we  know  that 


dx 


■p-2px  =  0  , 


(2) 


which  has  the  solution  x  =  x^  =  ll2.  Substituting  this  value  into  the  expres¬ 
sion  for  the  map  gives 


^max  =  /(^  =  ^p)  =  -  ^p)  =  4  '  (3) 

which  is  the  maximum  height  of  the  map.  Successive  iterations  of  the  map 
defined  by  equation  (1)  produce  a  sequence  of  points  x^,  which  wander 
ergodically  [4]  over  an  attracting  set  defined  by  a  minimum  and  maximum 
edge  (^max  defines  the  peak  of  the  map,  but  the  maxi¬ 

mum  boundary  of  the  attracting  set  as  well).  We  determine  the  minimum 
boimdary  of  the  attracting  set,  by 

^min  =  /(^  =  ^max)  =  "  ^max)  =  ^  ^^(4"  P)  ■  (4) 

This  value  is  important  in  determining  constraints  on  the  grammar  of  the 
map,  which  we  discuss  in  section  6.  Also  note  that  we  are  considering  only 
state  points  that  lie  in  the  unit  interval  on  the  real  axis.  This  is  true  for  0  <  p 
<4. 


(a) 


2x1  O'*  4x10'*  6x10'*  8x10'* 
Number  of  iterations,  n 


1x10^ 


Figure  2.  (a)  Return  map  and  (b)  sequential  output  (attracting  set)  of  logistic  map  for  p  =  3.97  in 
presence  of  uniformly  distributed,  zero-mean  noise  of  maximum  amplitude 
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3.  Symbolic  Dynamics  Representation  of  1-D  Chaotic 
Maps 


If  we  partition  the  unit  interval  into  disjoint  subsets  called  "atoms"  or 
"partition  elements"  so  that  a  state  point  falling  into  a  given  partition  ele¬ 
ment  is  assigned  a  particular  symbol,  then  we  form  a  symbolic  representa¬ 
tion  of  the  d)mamics.  For  a  binary  s5Tnbolic  system,  we  choose  the  peak  of 
the  map  (a:  =  1/2)  as  the  partition  boundary.  The  choice  of  this  partition  is 
not  arbitrary.  For  the  logistic  map,  this  point  is  known  as  the  generating  par¬ 
tition  because  it  and  all  its  iterates  and  pre-iterates  divide  the  state  space 
arbitrarily  finely  [16].  In  other  words,  each  state  point  has  its  own  unique 
symbol  sequence  assigned  to  it.  As  is  iterated  forward,  we  can  assign  a 
"0"  to  a  symbol  when  a  state-space  point  lands  on  the  left  side  of  the 
partition,  and  assign  a  "1"  to  the  symbol  when  a  state-space  point  lands 
on  the  right  side  (see  fig.  2a).  The  bits  are  shifted  into  the  right  side  of  an  N- 
bit  register  and  the  left-most  bit  is  thrown  out  (see  fig.  3).  For  example, 
a  16-bit  symbol  sequence  describing  a  trajectory  that  continuously  passes 
through  the  fixed  point  would  be  B„  =  1111111111111111  (the  fixed  point,  or 
period-1  point,  is  the  state  point  for  which  =  x„).  A  period-2  sequence 
could  be  =  1010101010101010  or  =  0101010101010101.  Now  we  want 
to  know  what  symbol  sequence  we  can  expect,  given  that  we  have  just 
passed  through  a  particular  point  in  state  space.  In  other  words,  we  ask 
what  state  point  sources  a  given  symbol  sequence.  In  order  to  do  this  we 
create  a  binary  coding  function  [6,7].  (From  this  point  on,  we  work  with  16- 
bit  symbol  sequences.)  We  define  the  number  k  as  the  value  of  the  symbol 
sequence  simply  by 


=  (5) 

Now  setting  N  =  16  and  iterating  the  map  16  times,  we  collect  the  set  of 
state  points  [xq,  Xj,  ...  Xjg]  and  the  corresponding  symbol  sequence  B  =  [feg, 
fcp  ...  fcjg].  We  say  that  the  first  state  point  Xg  evolves  the  symbol  sequence  B, 
which  has  value  k.  (In  the  following,  for  simplicity  of  notation,  we  refer  to 
k  interchangeably  as  the  value  of  the  symbol  sequence  and  the  sequence  it¬ 
self.)  Continuous  iteration  of  the  map  from  various  source  locations  will 
form  a  state-point /symbol-sequence  relationship  for  the  map  called  the  in¬ 
verse  coding  function: 

s[k]  =  Xj^ .  (6) 

That  is,  the  state  point  Xj^  sources  the  symbol  k. 

Forming  this  inverse  coding  function  is  part  of  what  we  call  the  character¬ 
ization  of  the  system.  Figure  4  shows  a  16-bit  inverse  coding  function  for 


Figure  3.  Bit-shifting 
diagram  for  a  finite- 
symbol  sequence 
register. 
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the  logistic  map  when  p  =  3.97.  Despite  the  imusual  topology  of  the  "func¬ 
tion,"  every  state  point  sources  a  unique  symbol  sequence.  We  will  see  that 
this  structure  is  based  primarily  upon  the  ordering  of  the  symbol  se¬ 
quences  across  the  state  space.  S)mibol  sequences  are  evolved  from  finite 
regions,  of  the  state  space  that  have  widths  on  the  order  of  2"^^.  *  A 
simple  way  of  determining  a  single  source  point  for  a  given  symbol  se¬ 
quence  k  is  by  averaging  over  the  state  points  producing  that  symbol 
sequence  (all  state  points  that  fall  into  a  subregion  will  produce  the  se¬ 
quence  k).  Namely, 


=  (7) 

where  N^.  is  the  number  of  observed  state  points  that  source  the  given  sym¬ 
bol  sequence  k  during  data  collection,  and  j  is  the/th  state  point  that  falls 
into  the  subregion  Aj^.  We  can  find  a  sequential  expression  for  the  partition 
refinements  if  we  first  realize  that  the  two-bit  partition  Xp2  occurs  where 
(one-bit  partition,  x^)  intersects  the  other  words,  we 

solve 


f(Xp2)  =  rXj,2  0--Xp2)  =  Xpi, 


(8) 


Figure  4. 16-bit  coding 
function,  s[k],  for  the 
logistic  map  at  p  = 

3.97  in  presence  of 
10“^^  noise. 
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*  Inverse  images,  that  is,  "backward"  iterations  of  the  partition  {[9, 1/2),  [1/2, 1)},  and  successive  intersections,  form 
refinements  of  the  partition  with  2N  pieces  of  average  width  =  2”^  at  the  Nth  stage  of  refinement. 
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which  is  a  quadratic  expression  from  which  we  can  derive  the  general  ex¬ 
pression  for  the  refinements,  * 


p±^p^-Apx^^ 
2p  ^ 


(9) 


If  there  is  a  complex  solution  to  equation  (9)  for  a  given  partition,  then  that 
finite-length  sequence  is  not  allowed;  thus  the  grammar  is  constrained. 
The  ordering  of  the  symbol  sequences,  across  the  map  from  left  to  right, 
conforms  to  a  Gray  ordering  [17].  Table  1  shows  three-bit  sequences  or¬ 
dered  according  to  the  standard  binary  "counting"  order,  and  ordered  ac¬ 
cording  to  a  "Gray  code." 

Figure  5  shows  a  plot  of  the  binary  value  of  sequences  in  a  Gray  ordering 
versus  their  position  in  the  list.  Note  the  similarity  to  the  coding  function 
in  figure  4.  If  the  partitions  are  known  and  the  ordering  is  known,  then  ev¬ 
erything  necessary  to  characterize  the  system  is  known.  In  other  words, 
given  the  expression  for  the  map,  a  complete  N-bit  coding  function  can  be 
calculated.  This  would  serve  to  significantly  reduce  the  characterization 
time.  For  systems  where  the  map  is  not  known  explicitly  (e.g.,  experimen¬ 
tal  systems,  complex  models),  yet  the  system  is  dissipative,  curve-fitting 
techniques  can  be  used  to  fit  the  return  map,  and  root-finding  methods  can 
be  used  to  determine  the  refinement  boundaries.  Figure  6  shows  a  four-bit 
refinement  completely  determined  by  this  method. 

Figure  7  shows  the  probability  mass  function  for  16-bit  symbol  sequences. 
It  is  apparent  that  certain  sequences  have  a  higher  probability  of  occur¬ 
rence  than  others.  There  are  also  constraints  on  the  system  that  disallow 
certain  sequences. 


Table  1.  Three-bit 
Gray  ordering 
example. 


k 

Binary 

Gray 

ordering 

0 

000 

000 

1 

001 

001 

2 

010 

on 

3 

oil 

010 

4 

100 

no 

5 

101 

111 

6 

110 

101 

7 

111 

100 

*The  partition  is  a  collection  of  intervals  P  =  IPil.  Refinement  is  P^  =  P  vf  ^  P  v..vf^  ^P  at  step  n,  where  P  vQ 
=  IPinQiWf  PiePandQieQ}. 
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Figure  6.  Four-bit 
partition  diagram  for 
the  logistic  map. 


Binary  sequence 


00000  0  00  11  1  111 
00001  1  11  11  1  100 
00111  1  00  00  1  111 


4.  Characterization  Under  the  Influence  of  Reference 
Perturbations 


We  now  deternvine  what  happens  to  the  system,  namely  the  inverse  cod¬ 
ing  function  (which  is  a  complete  description  of  the  symbolic  d)mamics), 
when  influenced  by  small  perturbations.  In  this  phase  of  the  characteriza¬ 
tion,  we  form  a  new  inverse  coding  function,  q[k],  which  describes  the  dy¬ 
namics  of  the  system  imder  the  influence  of  reference  perturbations.  As  the 
state  trajectory  generates  the  first  symbol,  we  add  a  small  reference  pertur¬ 
bation,  Pj.^,  to  the  parameter  p.  The  map  becomes 

=  [p  +  Pref] •  (10) 


The  next  15  symbols  are  observed  but  the  parameter  is  set  back  to  r.  We 
thus  see  the  effect  of  a  small  perturbation  on  the  evolving  sequence.  We 
now  have  a  new  inverse  coding  function  that  tells  us  what  state  point 
sources  each  allowed  sequence  if  a  reference  perturbation  is  present.  Figure  8 
shows  that  q[k]  is  qualitatively  similar  to  s[k],  shown  in  figure  4. 

The  reference  perturbation  actually  changes  the  symbolic /wfwre  of  the 
state-space  trajectory  (the  current  symbol  is  not  affected  as  the  perturba¬ 
tion  is  applied).  Consequently,  we  can  view  this  as  changing  the  effective 
sourcing  state  point  for  the  symbol  sequence  k.  As  a  result  of  our  charac¬ 
terization  phases,  we  know  that  for  a  reference  perturbation,  there  will 
be  a  measurable  difference  between  the  state  point  that  sources  the  symbol 
k  (i.e.,  q[k])  and  the  state  point  that  sources  the  symbol  k  under  zero  pertur¬ 
bation  (i.e.,  s[k]).  In  other  words. 


=  s[k]  -  q[k] . 


Figure  8. 16-bit  coding 
function,  q[k],  under 
influence  of  10“^ 
uniform  reference 
perturbations  applied 
every  16  iterations. 
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Using  this  result  we  form  a  proportionality  constant  for  each  symbol 
sequence; 


=  ^  -  (12) 

Figure  9  shows  a  plot  of  the  calculated  proportionality  constant  as  a  func¬ 
tion  of  the  symbol  sequence  k.  (Like  q,  d  has  fractal  structure.)  If  we  define 
the  initial  sequence,  or  source  sequence,  k^  =  k9  and  the  target  sequence,  kj  = 
k^^,  then  the  goal  is  to  go  from  k9  to  in  16  iterations.  We  make  incremen¬ 
tal  steps  by  bit-wise  shifting  k^^  into  ^k^  —>k^  •— >  k^^.  The  bits  from 

are  shifted  into  kP,  forming  each  subsequent  symbol  sequence.  Control 
perturbations  are  calculated  and  applied  so  that  a  given  state  point  sources 
a  desired  symbol  sequence,  rather  than  the  symbol  sequence  it  would  nor¬ 
mally  have  sourced  under  zero  perturbation.  Now  we  determine  the  ac¬ 
tual  state  point  that  has  occurred  via  iteration  of  the  map,  and  compare  it 
with  the  value  of  the  inverse  coding  function  that  will  source  the  desired 
symbol.  The  error  in  that  comparison  is  used  to  calculate  the  control 
perturbation: 


(13) 


Figure  9. 
Proportionality 
constant  d[k]. 
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5.  Targeting  the  Logistic  Map 

In  targeting,  we  are  trying  to  find  an  easy  route  from  one  state  point  to  an¬ 
other.  An  optimal  route  would  minimize  the  iterations  needed  to  get  from 
the  source  to  the  target,  as  well  as  minimize  the  control  energy  needed  to 
travel  there.  We  must  point  out  that  several  routes  will  be  possible  in  an 
ergodic  system.  If  we  wait  long  enough,  the  natural  dynamics  of  the  sys¬ 
tem  may  eventually  carry  a  source  point  near  a  target  point.  However, 
there  may  be  good  reasons  not  to  wait  that  long  [11].  Targeting  by  control¬ 
ling  symbolic  d5mamics  provides  a  navigational  method.  The  number  of  it¬ 
erations  needed  to  get  to  within  a  neighborhood  of  the  target  and  the  pertur¬ 
bations  needed  to  guide  the  trajectories  there  are  completely  determined 
by  the  number  of  bits  in  the  symbol  sequence  register.  Once  we  realize  that 
a  point  in  state  space  is  uniquely  linked  to  a  symbol  sequence,  targeting 
becomes  a  simple  matter  of  shifting  in  a  target  sequence  (kj)  from  a  source 
sequence  (k^).  The  procedure  is  as  follows: 

•  Iterate  the  system  without  and  with  perturbation  to  obtain  s[A:]  and  q[k]. 

•  Define  a  source  point  x^. 

•  Search  s[k]  in  order  to  find  the  closest  value  to  s[/t5]  =  x^. 

•  Define  a  target  point  Xj. 

•  Search  s[k]  in  order  to  find  the  closest  value  to  xf.  slfcj]  =  Xj. 

•  Sequentially  shift  in  the  bits  from  kj  into  kg.  Calculate  and  apply  the  neces¬ 
sary  perturbations  for  each  successive  shift  by  equation  (13). 

The  first  example  is  targeting  an  arbitrary  state  point  from  an  arbitrary 
source  point  Xg-^Xj  =  Xg-^ Xj.  In  this  example,  we  demonstrate  control 
of  a  period-16  orbit  by  simply  targeting.  Table  2  shows  the  results. 

Source  symbol  kg  =  28574  =  D[0110111110011110]  and  target  point  kj  = 
56466  =  D[1101110010010010],  where  D[  ]  takes  the  decimal  value  of  the 
quantity  in  the  square  brackets.  Figure  10  shows  the  path  taken  on  the 
attractor  graphically. 

A  second  example  is  what  we  would  term  a  rapid  period-1  stabilization. 
From  an  arbitrary  source  point,  we  target  the  fixed  point  {kj  =  65535  = 
D[llllllllllllllll])  and  then  repeatedly  target  the  fixed  point.  Table  3 
shows  the  results  of  this  example. 

In  both  examples,  the  control  was  carried  out  in  the  presence  of  uniformly 
distributed  noise  of  10”^^  peak-to-peak  amplitude  and  zero  mean.  The 
magnitude  of  the  noise  has  a  large  effect  upon  the  size  of  the  perturbations 
needed  to  target  and  stabilize  the  system.  Since  the  range  of  the  map  is  ap¬ 
proximately  unity,  the  size  of  the  perturbations  shown  in  table  3  is  good 
for  general  comparison.  Figure  11  shows  the  results  graphically. 
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Table  2.  Arbitrary  state-point  targeting  and  demonstration  of  a  period-16  control  sequence. 


Iteration 

i 

Current 

symbol 

k 

Associated 
state  point* 
s[fc] 

Actual  Difference  between 

state  point,  state  and  target  points, 
x{n+l)  \x{n+l)-Xj\ 

Control 

bit, 

bi 

Applied 

perturbation. 

Pi 

Source,  Xr 

28574 

1.966511E-01 

1.966377E-01 

4.388792E-01 

— 

— 

0  ' 

57149 

6.2715820341E-O1 

6.2714610232E-01 

8.3707726855E-03 

1 

-5.2863171234E-05 

1 

48763 

9.2830466935E-01 

9.2830809745E-01 

2,9279122244E-4)1 

1 

1.7791012268E-04 

2 

31990 

2.6425563888E-01 

2.6422396978E-4)1 

3.7129290523E-O1 

0 

3.0456946497E-04 

3 

63981 

7.7184712831E-01 

7.7186557569E-01 

1.3634870068E-01 

1 

2.2956528329E-04 

4 

62427 

6.9910684008E-01 

6.9911418577E-01 

6.3597310759E-02 

1 

5.5389791500E-05 

5 

59319 

8.3509921281E-01 

8.3511520937E-4)1 

1.9959833436E-01 

1 

3,1054768026E-04 

6 

53102 

5.4673072082E-01 

5.4670301369E-01 

8.8813861317E~02 

0 

^.4358120314E-4)5 

7 

40668 

9.8382616628E-01 

9.8382975647E-01 

3.4831288146E-01 

0 

8,7426558587E-04 

8 

15801 

6.3167764857E-02 

6.3171712545E-02 

5.7234516247E-4)1 

1 

-2.3937827438E-04 

9 

31602 

2.3492524516E-01 

2.34934591  lOE^Il 

4.0058228391E-01 

0 

-1.0963326411E-f)4 

10 

63204 

7.1352847336E-01 

7.1354940069E-O1 

7.8032525677E-02 

0 

1.7461076372E-04 

11 

60873 

8.1148216065E-01 

8.1149040418E^1 

1.7597352917E--01 

1 

1.3262176675E-04 

12 

56210 

6.0733884086E-01 

6.0732598829E-01 

2.8190886722E-02 

0 

-4.6614076940E-05 

13 

46884 

9.4674932521E^1 

9.4675897847E-01 

3.1124210346E-01 

0 

6.8965761921E-04 

14 

28233 

2.0014128684E-01 

2.0014823141E-01 

4.3536864360E-O1 

1 

-1.0344111909E-O4 

15 

56466 

6.3551220592E-01 

6.3553644028E-01 

1.9565267675E-05 

0 

1.1374004902E-04 

ixr  =  xs) 

16 

47397 

9.1958692689E-01 

9.1959694288E~01 

2.8408006787E^1 

1 

4.4703382471E^4 

17 

29259 

2.9354324144E-01 

2.9356852297E-01 

3.4194835204E-01 

1 

-2.0087323557E-04 

18 

58518 

8.2326609079E-01 

8.2328094156E^1 

1.8776406655E^1 

0 

2.6642402486E-04 

19 

51501 

5.7761373061E-01 

5.7763181026E^1 

5.7885064750E-02 

1 

4.6868864940E-05 

20 

37467 

9.6857807135E-01 

9.6858544389E-01 

3.3306856888E-01 

1 

9.0439706514E-04 

21 

9399 

1.2081458035E-01 

1.2082541540E-01 

5.1469145961E-01 

1 

-3.0801 158981E-04 

22 

18798 

4.2166424384E^1 

4.2168701957E-4)1 

2.1382985544E~01 

0 

-5.9718902906E-05 

23 

37596 

9.6813298861E-01 

9.6813773265E-01 

3.3262085764E-01 

0 

5.7208161307E-04 

24 

9657 

1.2247435751E--01 

1.2248048826E*-01 

5.1303638675E-4)1 

1 

-1.6976738140E-04 

25 

19314 

4.2664417153E-01 

4.2667345610E^1 

2.0884341891E-01 

0 

-7.1813376103E^5 

26 

38628 

9.7113015096E-01 

9.7113660814E-01 

3.3561973313E-01 

0 

8.6391995156E-04 

27 

11721 

1.1129923533E-01 

1.1130449298E-01 

5.2421238203E-01 

1 

-1.6571221563E-04 

28 

23442 

3.9265628538E-01 

3.9267934570E-01 

2.4283752931E-01 

0 

-8.3757283727E-05 

29 

46884 

9.4674932521E-01 

9.4675466576E--01 

3.1123779075E-01 

0 

3.8154485518E-04 

30 

28233 

2.0014128684E-01 

2.0014800027E-01 

4.3536887474E-01 

1 

-9.9998254070E-05 

31 

56466 

6.3551220592E-01 

6.3553644116E-01 

1.9566145705E-05 

0 

1.1374416990E-04 

*  Stored  in  unperturbed  coding  function. 


Figure  10.  Graphical 
description  of 
(a)  targeting  an 
arbitrary  state  point 
from  an  arbitrary 
source  state  point, 
and  (b)  a  period-16 
orbit 
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Table  3.  An  example  of  rapid  fixed-point  stabilization. 


Iteration  Current  Associated  Actual  Difference  between  Control  Applied 

i  symbol  state  point*  state  point,  state  and  target  points,  bit,  perturbation, 

k  s[/c]  x{n+l)  \  x{n-\-l)  - Xj\ 


Source,  x^ 

24524 

0 

49049 

1 

32563 

2 

65127 

3 

64719 

4 

63903 

5 

62271 

6 

59007 

7 

52479 

8 

39423 

9 

13311 

10 

26623 

11 

53247 

12 

40959 

13 

16383 

14 

32767 

15  =  Xr 

65535 

16 

65535 

17 

65535 

18 

65535 

19 

65535 

31 

65535 

3.754551E-01 

9.3092592112E-01 

2.5531107999E-01 

7.5477741050E-01 

7.3478575310E-01 

7.7363884542E-01 

6.9520811550E~01 

8.4119409510E-01 

5.3031158533E-01 

9.8884884466E-01 

4.3769333556E-02 

1.6617644513E-01 

5.5006584947E-01 

9.8254512851E-01 

6.8080295032E-02 

2.5189261197E^)1 

7.4810942348E-01 

7.4810942348E-01 

7.4810942348E-01 

7.4810942348E-O1 

7.4810942348E-01 

7.4810942348E-01 


3.754377E-01 

9.3090239932E-01 

2.5528002008E-01 

7.5480580070E-01 

7.3480194841E-01 

Z7365656806E-01 

6.9523387801E-4)1 

8.4121846190E-01 

5.3033912777E-01 

9.8885339495E-01 

4.3776756851E-02 

1.6615818535E-01 

5.5009060092E-01 

9.8254944157E-01 

6.8086293661E-02 

2,5187818663E~01 

7.481 1771 656E-01 

7.4811367768E-01 

7.4811364222E-01 

7.4811364191E^1 

7.4811364191E-01 

7.481136419iE-01 


3.726731E-01  - 

1.8279157431E-^1  1 

4.9283080493E>-01  1 

6.6949756866E-03  1 

1.3308876602E-O2  1 

2.5545743052E-02  1 

5.2876946997E-4)2  1 

9.3107636892E-02  1 

2.1777169724E-01  1 

2.4074256994E-01  1 

7.0433406816E-01  1 

5.8195263966E-01  1 

1.9802022409E-01  1 

2.3443861656E-01  1 

6.8002453135E-01  1 

4.9623263838E-01  1 

6.8915503201E^6  1 

2.8526694139E-06  1 

2.8172131984E-06  1 

2.8169006118E-06  1 

2.8168978559E-06  1 

2.8168978313E-06  1 


-1.2868674989E-03 

3.1896364866E-04 

3.1421345006E-04 

1.5615314549E-04 

2.2219717700E-04 

1.8898253249E-04 

4.9697939082E-04 

3.0639647084E-05 

1.6325515489E-03 

-6.5488648724E^4 

3.5019205078E-04 

4.2194143223E-05 

9.6400954124E-O4 

-3.2305512045E-04 

1.5135440551E-04 

8.7088069641E-4)5 

4.4674574484E-05 

4.4302238171E-05 

4.4298955605E-05 

4.4298926665E-05 

4.4298926407E-05 


Figure  11.  Graphical 
description  of  rapid 
fixed-point 
stabilization. 
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6,  Constrained  Sequences 

If  we  consider  targeting  as  equivalent  to  navigating  routes  from  an  origin 
to  a  destination,  then  symbolic  dynamics  provides  a  formalized  road  map 
for  the  journey.  As  is  evident  by  figure  7,  there  are  several  symbol  se¬ 
quences  that  have  either  a  very  low  or  zero  probability  of  occurrence. 
These  constraints  in  the  symbol  space  arise  because  the  map,  for  0  <  r  <  4, 
does  not  extend  over  the  whole  range  [0, 1]  (see  fig.  6).  In  fact,  from  equa¬ 
tion  (4)  we  determine  that  =  0.029551687512,  and  thus  the  minimal  se¬ 
quence  derived  by  iterating  the  map  forward  and  collecting  the  s)mibols  is 
Bmin  =  0001011001001011.  The  minimal  sequence  is  the  smallest  sequence 
in  magnitude,  in  a  Gray  counting  order,  possible  by  the  map.  Table  4  dem¬ 
onstrates  this  iteration. 


Had  fallen  exactly  on  one  of  the  refinement  botmdaries  (see  fig.  6),  we 

would  have  had  a  simple  zero  run-length  constraint  (for  a  run-length  con¬ 
straint  on  zeroes  alone,  no  more  than  K  zeroes  can  occur  in  a  row).  The 
state  points  would  have  mapped  to  peak  in  K  iterations.  Our  constraint  for 
r  =  3.97  is  much  more  complicated  than  this.  However,  by  viewing  the 
minimal  sequence,  we  do  know  that  no  symbol  sequence  can  occur  that 
has  more  than  three  zeroes  in  a  row  within  it. 


A  measure  of  the  chaos  of  the  map  can  be  quantified  in  a  measure  called 
the  topological  entropy  [18],  hj.  In  general,  this  entropy  is  an  information- 
theoretical  quantity  that  provides  a  measure  of  the  possible  states  that  a 
system  can  take  on.  In  the  context  of  chaotic  dynamics,  we  can  use  this 
measure  to  quantify  the  number  of  orbits  possible  in  a  chaotic  system. 
Since  orbits  in  state  space  are  linked  to  symbol  sequences  in  this  symbolic 
dynamics  formalism,  we  can  approximate  this  entropy  by  comparing  the 
number  of  N-bit  symbol  sequences  generated  by  the  system  versus  all  pos¬ 
sible  N-bit  symbol  sequences.  This  measure  has  a  positive  value  for  chaotic 


Table  4.  Iteration  of 
^min  yielding  minimal 
sequence  = 
0001011001001011 
(5707). 


Iteration 

n 

Point 

Corresponding 

symbol, 

K 

0 

0.0295516875  12 

0 

1 

0.113853189575 

0 

2 

0.400535843758 

0 

3 

0.953224320043 

1 

4 

0.177013231413 

0 

5 

0.578347802852 

1 

6 

0.968130638499 

1 

7 

0.122489210033 

0 

8 

0.426717845731 

0 

9 

0.971180011687 

1 

10 

0.111117904449 

0 

11 

0.392119741566 

0 

12 

0.946296543866 

1 

13 

0.201752997885 

0 

14 

0.639363441146 

1 

15 

0.915393990149 

1 
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systems  and  a  value  of  zero  for  nonchaotic  systems  [4].  For  the  character¬ 
ization  that  we  carried  out  with  one  million  iterations  of  the  map,  approxi¬ 
mately  40  percent  of  the  possible  65,536  sequences  were  allowed.  The 
measure  hj  is  only  approximated  for  a  finite-length  symbolic  dynamics 
representation.  It  is  calculated  as 


hr=  lim 


log2n  _  log2n 
N  “  16 


=  0.92  , 


(14) 


where  n  =  the  number  of  allowed  sequences  and  N  =  the  number  of  bits. 
Other  disallowed  sequences  were  due  to  more  complex  symbol  constraints 
and  possibly  occurred  because  of  sequences  that  have  a  low  probability  of 
occurrence. 

The  constraint  on  symbol  sequences  presents  two  problems:  (1)  some  re¬ 
gions  may  not  be  able  to  be  targeted  simply  because  they  are  not  on  the 
attractor  and,  most  importantly,  (2)  derived  routes  may  not  be  passable 
because  they  involve  shifting  in  disallowed  sequences.  The  second  case 
was  not  at  all  uncommon.  However,  because  the  attractor  exists,  any  exist¬ 
ing  state  point  on  the  attractor  is  accessible  to  any  other  existing  state 
point.  The  route  may  be  longer,  but  a  route  does  exist. 

One  method  of  dealing  with  problem  2  is  the  insertion  of  a  buffer  sequence 
of  length  Ng  between  the  source  sequence  and  the  target  sequence  (see  fig. 
12).  The  buffer  sequence  may  or  may  not  be  as  long  as  the  source  and  tar¬ 
get  sequences.  For  example,  if  there  is  only  a  simple  run-length  constraint 
(i.e.,  any  sequence  containing,  say,  four  zeroes  or  more  in  a  row  is  disal¬ 
lowed),  then  inserting  the  one-bit  buffer-sequence  "1"  will  ensure  that 
runs  of  more  than  three  zeroes  will  not  be  generated  by  concatenation.  If 
Ng  =  N,  we  call  this  buffer  sequence  a  jump  point,  because  it  directly  trans¬ 
lates  to  a  point  in  state  space.  It  may  be  possible,  for  some  systems,  to  find 
a  single  jump  point  in  state  space  that  will  act  as  a  conduit  to  and  from  all 
possible  state  points  (i.e.,  the  fixed  point).  If  a  single  jump  point  does  not 
exist,  it  may  be  possible  to  create  a  table  of  buffer  words  that  will  allow 
safe  passage  from  varied  source/target  combinations.  Table  5  shows  an 
example  of  an  attempt  to  target  via  an  unpassable  route  (road  block),  and 
table  6  shows  the  insertion  of  a  jump-point  sequence  to  allow  targeting  to 
occur.  Instead  of  16  iterations,  there  are  32.  The  following  bit  shifting  dia¬ 
gram  shows  a  configuration  in  which  a  disallowed  sequence  of  four  zeroes 
would  occur. 

1111111111111100  ^  0011101110010010 

Source  sequence  Target 


Figure  12.  Sequence  -  < - 

shifting  diagram 

showing  insertion  of  I  I  | - 1  - 

a  buffer  word.  I _ I  | _ |  _ 

ks  kB  kj 
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Table  5.  Attempt  to  target  via  an  unpassable  route  (roadblock):  target  =  7.627758E-02. 


Iteration 

i 

Current 

S5nTibol 

k 

Associated 
state  point* 
s[k] 

Actual  Difference  between  Control  Applied 

state  point,  state  and  target  points,  bit,  perturbation, 

x{n+l)  \x{n+l)-Xj\ 

Source,  Xc 

:  65532 

7.480572E-01 

7.48572E-01 

6.1718012E-01 

— 

— 

0 

65528 

7.4823513083E-01 

7.4821648997E-4)1 

67193890867E-01 

0 

-1.9637686395E-04 

1 

65520 

-l.OOOOOOOOOOE-03 

7.4786570304E-01 

6.7158812173E-01 

0 

O.OOOOOOOOOOE+00 

2 

65505 

-1.0000000000E-<)3 

7.4859355423E-01 

6.7231597292E-01 

1 

O.OOOOOOOOOOE+00 

3 

65475 

-l.OOOOOOOOOOE-03 

7.4715900084E--01 

6.7088141953E-01 

1 

O.OOOOOOOOOOE+00 

4 

65415 

-l.OOOOOOOOOOE-03 

7.4998239936E-01 

6.7370481806E-01 

1 

O.OOOOOOOOOOE+00 

5 

65294 

-1.0000000000E~03 

7.4440999502E-01 

6.6813241372E-01 

0 

O.OOOOOOOOOOE+00 

6 

65053 

-l.OOOOOOOOOOE-03 

7.5534716369E-*01 

6,7906958239E-01 

1 

o,ooooooooooE+oo 

7 

64571 

-l.OOOOOOOOOOE-03 

7.3364742818E-01 

6.5736984688E-4)1 

1 

O.OOOOOOOOOOE+00 

8 

63607 

-l.OOOOOOOOOOE-03 

7.7577330983E-01 

6.9949572853E--01 

1 

O.OOOOOOOOOOE+00 

9 

61678 

-l.OOOOOOOOOOE-03 

6.9057791289E-01 

6,1430033159E--01 

0 

O.OOOOOOOOOOE+00 

10 

57820 

-l.OOOOOOOOOOE-03 

8.4830989369E>Ol 

7.7203231239E-431 

0 

O.OOOOOOOOOOE+00 

11 

50105 

-l.OOOOOOOOOOE-03 

5.1086052428E-01 

4.3458294298E-01 

1 

O.OOOOOOOOOOE+00 

12 

34674 

-l.OOOOOOOOOOE-433 

9.9203179357E-01 

9.1575421227E-01 

0 

O.OOOOOOOOOOE+00 

13 

3812 

-l.OOOOOOOOOOE-03 

3.1381774037E-02 

4.4895807264E-02 

0 

O.OOOOOOOOOOE+00 

14 

7625 

-1.0000000000E^3 

1.2067598342E-01 

4.4398402120E-02 

1 

O.OOOOOOOOOOE+OO 

15 

15250 

7.6277581300E-02 

4.2126982206E-01 

3.4499224076E-01 

0 

-1.6979509865E+01 

20 


Table  6.  Insertion  of  a  jump-point  sequence:  target  =  7.627758E-02. 


Iteration 

Current 

Associated 

Actual 

Difference  between 

Control 

Applied 

i 

symbol 

k 

state  point* 
s[fc] 

state  point,  state  and  target  points, 
x{n+l)  \  x(n-\-l)  -  Xj\ 

bit, 

bi 

perturbation. 

Pi 

Source, 

65532 

7.480572E-01 

7.480572E-01 

6.7178012E-01 

— 

— 

0 

65529 

7.4821534803E-01 

7.4821648997E-01 

6.7193890867E-01 

1 

1.2020856267E-05 

1 

65523 

7.4789678410E-01 

7,4790496279E-01 

6.7162738148E-4)1 

1 

8.5705721635E-05 

2 

65511 

7.4851148670E-01 

7.4853244204E^1 

6.7225486074E-01 

1 

2.1938233599E-4)4 

3 

65487 

7.4730639913E-4)1 

7.4732090595E-01 

6.7104332464E-01 

1 

1.5194295363E-04 

4 

65439 

7.4968037285E-01 

7.4969325758E-01 

6.7341567628E-01 

1 

l,3654124576E-04 

5 

65343 

7.4498850371E-4)1 

7.4500919156E-01 

6.6873161026E-01 

1 

2.1259724099E--04 

6 

65151 

7.5420408066E-01 

7.5422331552E-01 

6.7794573421E-01 

1 

2.1137871668E-4)4 

7 

64767 

73594964817E-01 

7,3596015059E-4)1 

6,5968256928E^1 

1 

1.0186318458E-04 

8 

63999 

7.7146471177E-01 

7.7148139845E-01 

6.9520381715E-01 

1 

2.0284115455E-O4 

9 

62463 

6.9991503755E-01 

6.9993828512E-01 

6.2366070382E-01 

1 

1.7779241532E-04 

10 

59391 

8.3382614509E-01 

8.3383538790E-01 

7.5755780659E-01 

1 

1.7822670304E-04 

11 

53247 

5.5006584947E-01 

5.5008387014E-01 

4.7380628884E-01 

1 

3.0720087778E-05 

12 

40959 

9.8254512851E-01 

93254933758E-01 

9.0627175628E-01 

1 

9.4076647475E-04 

13 

16383 

6.8080295032E-02 

6.8086352596E-02 

8.1912287050E-03 

1 

-3.2622907802E-04 

14 

32767 

2.5189261197E~01 

2,5187824627E-01 

1.7560066497E-01 

1 

1.5072867990E--04 

15 

65535 

7.4810942348E-01 

7.4811777508E-01 

6.7184019377E-01 

1 

8.7702556859E-05 

(jump  point) 
16 

65534 

7.4812887863E-01 

7.4811373712E-01 

6.7183615582E-01 

0 

-1.5891211003E-^ 

17 

65532 

7.4805721716E-01 

7.4807522026E-01 

6.7179763896E-01 

0 

1.8914202384E-4)4 

18 

65529 

7.4821534803E-01 

7.4821668426E-01 

6.7193910296E-01 

1 

1.4066086278E-05 

19 

65523 

7.4789678410E-01 

7.4790496517E-01 

6.7162738387E-01 

1 

8.5730697460E-05 

20 

65511 

7.4851148670E-01 

7.4853244205E-01 

6.7225486075E-01 

1 

2.1938251974E-4)4 

21 

65486 

7.4733673283E-01 

7.4732090595E-01 

6.7104332464E-01 

0 

-1.6619514188E-04 

22 

65437 

7.4964850806E^1 

7.4963318281E-01 

6.7335560151E-4)1 

1 

-1.6253677605E-04 

23 

65339 

7.4508423833E-01 

7.4507215149E-01 

6.6879457019E-01 

1 

-1.2437608627E-04 

24 

65143 

7,5405302280E-~01 

7,5403680808E-01 

6.7775922678E-01 

1 

-1.7810303433E-f)4 

25 

64750 

7.3625517555E~01 

7.3626426862E-01 

6.5998668732E-01 

0 

8.7520937169E-4)5 

26 

63964 

7.7091680690E-01 

7.7090845937E-01 

6.9463087807E^1 

0 

-1.0179477427E-04 

27 

62393 

7.0112712755E-01 

7.0111824952E-01 

6.2484066822E-01 

1 

--6.8284691909E-05 

28 

59250 

8.3189343770E-01 

8.3190500518E-01 

7.5562742388E-01 

0 

2.1789570484E-04 

29 

52964 

5.5521428967E-01 

5.5519162743E-01 

4.7891404613E-01 

0 

-A225SS1S703E-05 

30 

40393 

9.8039905049E-01 

9.8039654353E-4)1 

9.041 1896223E-01 

1 

-4.9230768205E-04 

31 

i^j) 

15250 

7.6277581300E-02 

7.6290666287E-02 

1.3084986331E-^5 

0 

-6.4400482124E-04 
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7.  Two-Dimensional  Maps 

More  general  three-dimensional  systems  can  yield  more  complicated  two- 
dimensional  Poincare  maps  and  thus  two-dimensional  return  maps.  An 
example  of  a  two-dimensional  map  is  generated  by  iteration  of  the 
expression 


(15) 


This  map  is  called  the  Henon  map  [4].  Figure  13  shows  a  plot  of  the  strange 
attractor  (or  fractal  attracting  set)  that  results  when  A  =  1.4  and  B  =  0.3. 
Proper  characterization  of  this  system  for  symbolic  dynamics  targeting  de¬ 
pends  upon  determining  the  correct  generating  partition.  It  has  been 
shown  that  the  proper  generating  partition  exists  where  the  Kolmogorov  en¬ 
tropy  (or  metric  entropy,  h^)  is  maximized  [5].  This  is  a  particularly  diffi¬ 
cult  measure  to  calculate,  especially  with  experimental  data.  In  special 
cases  where  the  metric  entropy  and  the  topological  entropy  are  equal,  find¬ 
ing  the  pomt(s)  where  hj  is  maximized  will  produce  the  generating  parti¬ 
tion.  Figure  14  shows  the  result  of  an  approximation  of  the  generating  par¬ 
tition  for  the  Henon  map  where  this  approach  was  used  with  16-bit 
S)anbolic  sequences. 


Figure  13.  Strange  attractor  for  two- 
dimensional  Henon  map  (A  =  1.4,  B  =  0.3). 


Topological  entropy,  hj 

Figure  14.  Topological  entropy  (/ij-) 
approximation  versus  horizontal  partition 
position  y  (showing  peak  occurring  at  y^  = 
-0.086). 
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8.  Control  Through  Crisis:  Ternary  Systems 

The  logistic  map  undergoes  a  crisis  when  p  >  4.  In  this  case  a  region  de¬ 
fined  by 


Figure  15.  Logistic 
map  undergoing 
crisis  for  p  >  4. 


will  project  the  map  to  a  point  greater  than  1,  and  consequently  send  it  off 
to  -oo  (see  fig.  15).  If  it  is  possible  to  characterize  the  system  near,  but  not  at, 
crisis,  then  it  is  possible  to  form  coding  functions  based  upon  a  ternary 
s5mibol  structure.  For  example,  "0"  is  assigned  to  any  point  passing  to  the 
left  of,  "1"  to  any  point  passing  to  the  right  of,  and  "2"  to  any  point  passing 
through  the  region  (this  method  would  demand  that  the  maximum  crisis 
region  be  known  a  priori).  Simply  stated,  any  region  of  state  space  that 
leads  to  crisis  can  be  given  the  sjmibol  "2."  If  the  system  is  targeted,  or  con¬ 
trolled,  through  only  the  sequences  containing  O's  and  I's  (thus  binary), 
then  even  when  p  >  4,  the  trajectories  will  never  enter  the  crisis  region. 

For  control  to  be  maintained  during  crisis,  the  near-crisis  characterization 
must  be  extrapolated.  This  brings  to  bear  the  notion  of  controlling  a  chaotic 
system  xmder  d)mamic  parameter  conditions.  For  a  system  where  the  de¬ 
scription  of  the  dynamics  is  known  explicitly,  such  as  the  logistic  map,  it 
would  be  possible  to  project  the  variations  of  the  symbolic  sequences  due 
to  varying  p.  Unperturbed  and  perturbed  coding  functions  can  be  deter¬ 
mined  for  various  values  of  p.  For  systems  where  there  is  no  explicit  ana¬ 
lytical  definition  of  the  map,  it  may  be  possible  to  curve-fit  and  param¬ 
etrize  the  map.  Rapid  characterization  could  then  be  calculated  by  the 
method  described  in  section  3.  A  single-humped  map  with  a  positive  first 
derivative  in  the  "0"  partition  will  have  a  Gray  ordering.  Regions  other 
than  those  that  induce  crisis  can  be  disallowed  by  judicious  choice  of  the 
ternary  partitions. 
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9.  Conclusion 


We  have  shown  that  the  formalism  developed  to  control  symbolic  dynam¬ 
ics  in  chaotic  systems  using  small  perturbations  readily  lends  itself  to  tar¬ 
geting.  An  important  application  that  could  be  affected  by  this  characteris¬ 
tic  is  rapid  stabilization  of  collapsible  systems.  Studies  of  voltage  collapse 
models  for  power-delivery  networks  have  shown  that  chaotic  oscillations 
often  precede  collapse  [11].  Systems  of  this  type  are  undergoing  a  form  of 
crisis.  An  algorithm  of  this  type  can  be  applied  to  quickly  guide  the  trajec¬ 
tories  of  these  systems  to  regions  of  state  space  where  periodic  stabiliza¬ 
tion,  using  the  OGY  principle  [3],  can  be  performed. 

This  targeting  method  readily  lends  itself  to  digital  processing  techniques. 
The  coding  functions  themselves  are  arrays  whose  storage  space  depends 
largely  upon  the  number  of  bits  defined  for  the  symbol  sequences.  Once 
characterization  is  complete,  only  s[k]  and  d[k]  need  be  stored,  which  are 
double-precision  arrays  2^  long,  where  N  is  the  number  of  bits  in  a  symbol 
sequence.  One  could  even  quantize  the  state  points  in  order  to  deal  strictly 
with  integer  quantities.  This  further  reduces  needed  storage  space  and 
speeds  up  calculation  time.  For  each  iteration,  the  control  perturbations 
are  calculated  by  a  subtraction  and  multiplication  (see  eq  (13)).  Given 
the  proper  generating  partition,  characterization  and  targeting  for  two- 
dimensional  maps  will  proceed  in  the  same  manner  as  for  a  one¬ 
dimensional  map. 

Although  we  know  the  analytical  expression  for  the  logistic  map,  it  was 
not  necessary  to  know  it  in  order  to  target  the  system.  A  long  enough  pe¬ 
riod  of  observation  and  learning  is  sufficient  to  create  a  model  from  measure¬ 
ment.  This  is  the  approach  that  many  researchers  and  experimenters  con¬ 
tinue  to  have  success  with. 
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